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LONG-TERM  GOALS 

1,  Passive  Source  Localization  in  a  Fluctuating  Ocean  Waveguide 

The  long-term  goal  of  this  research  is  to  determine  necessary  eonditions  on  eontrollable  experimental 
quantities,  such  as  sample  size  or  signal-to-noise  ratio  (SNR),  so  that  source  loealization  estimates  in  a 
fluetuating  ocean  waveguide  ean  attain  design  error  thresholds.  Using  asymptotie  expansions  of  the 
bias  and  covariance  of  maximum  likelihood  estimates  (MLEs)  in  inverse  orders  of  sample  size  or 
SNR,^  we  aim  to  derive  an  analytical  expression  for  the  second-order  covariance  of  MLEs  obtained 
from  general  eomplex  Gaussian  data  vectors.  This  expression  can  then  be  used  to  speeify  necessary 
conditions  for  accurate  estimation  in  many  practical  problems  where  both  the  mean  and  varianee  of  the 
measurement  are  functions  of  the  estimation  parameters.  Currently  sueh  eonditions  ean  only  be 
determined  in  eases  where  either  the  mean  or  the  variance  of  the  measurement  is  parameter  dependent, 
but  not  both.^  Our  ultimate  goal  is  to  quantify  the  effects  of  internal  waves  on  the  biases  and  varianees 
of  source  localization  estimates,  as  well  as  on  the  sample  sizes  and  SNRs  required  to  obtain  accurate 
estimates. 

2,  Velocity  Estimation  of  Underwater  Target  Swarms 

The  long-term  goal  of  this  researeh  is  to  develop  a  method  for  instantaneously  estimating  the  means 
and  standard  deviations  of  the  velocity  and  position  of  groups  of  self-propelled  underwater  seatterers 
observed  by  a  long-range  remote  sensing  system.  We  aim  to  obtain  these  estimates  based  on  an 
analytical  expression  for  the  magnitude  squared  of  the  range-velocity  ambiguity  function  of  the 
acoustic  field  seattered  from  such  target  groups.  This  expression  ean  then  be  used  to  estimate  the  first 
and  seeond  moments  of  the  targets’  veloeity  and  position.  Currently  no  technique  exists  to 
simultaneously  estimate  the  seeond-order  veloeity  and  position  statistics  of  underwater  target  groups. 
Our  ultimate  goal  is  to  utilize  such  estimates  to  instantaneously  distinguish  and  classify  underwater 
targets. 


1 


Report  Documentation  Page 

Form  Approved 

0MB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
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OBJECTIVES 


The  primary  objectives  of  the  research  are  to: 

•  Derive  an  analytical  expression  for  the  second-order  covariance  of  MLEs  obtained  from  general 
complex  Gaussian  data  vectors  when  both  the  measurement  mean  and  variance  are  parameter 
dependent. 

•  Determine  conditions  on  sample  size  and  SNR  so  that  source  localization  estimates  attain  desired 
error  thresholds  in  a  fluctuating  ocean  waveguide. 

•  Quantify  the  effects  of  ocean  fluctuations  on  the  biases  and  errors  of  parameter  estimates. 

•  Derive  an  analytical  expression  for  the  range-velocity  ambiguity  function  of  the  acoustic  field 
scattered  from  underwater  target  groups. 

•  Develop  a  method  for  distinguishing  and  classifying  groups  of  underwater  targets  by  estimating  their 
second-order  velocity  and  position  statistics. 

•  Publish  work  already  completed  on  optimally  estimating  Lambertian  surface  orientation  with 
application  to  high  frequency  sonar  imaging. 

APPROACH 

To  determine  the  conditions  necessary  to  obtain  accurate  source  localization  estimates  in  a  fluctuating 
ocean  waveguide,  we  first  expand  the  bias  and  covariance  of  MLEs  in  inverse  orders  of  sample  size  or 
SNR,  where  the  first-order  term  in  the  covariance  expansion  is  the  Cramer-Rao  lower  bound  (CRLB).^ 
We  then  derive  an  analytical  expression  for  the  second-order  covariance  of  MLEs  obtained  from 
general  complex  Gaussian  data  vectors.  Necessary  conditions  on  sample  size  and  SNR  are  determined 
by  requiring  that  the  MLEs  of  source  range  and  depth  (i)  become  asymptotically  unbiased  and  have  a 
variance  that  attains  the  CRLB,  and  (ii)  have  a  CRLB  that  is  within  the  desired  error  thresholds.  The 
first  condition  is  achieved  when  the  first-order  bias  and  the  second-order  covariance  are  smaller  than 
the  true  value  of  the  parameter  and  the  first-order  covariance,  respectively.'  To  model  the  waveguide 
containing  the  random  internal  waves  we  employ  an  already  available  physical  model  that  provides 
expressions  for  the  mean,  mutual  intensity,  and  spatial  covariance  of  the  acoustic  field  forward 
propagated  through  random  3-D  internal  waves  in  a  stratified  ocean  waveguide  for  a  continuous  wave 
narrowband  signal.^ 

To  estimate  the  instantaneous  velocity  and  position  statistics  of  underwater  target  groups,  we  derive 
analytical  expression  for  the  magnitude  squared  of  the  range-velocity  ambiguity  function  in  terms  of 
the  targets’  velocity  and  position  statistical  distributions.  This  expression  is  based  on  an  analytical 
model  for  the  acoustic  field  scattered  from  a  simple  harmonic  source  by  a  single  moving  target  in  a 
stratified  waveguide."'  We  then  show  that  the  moments  of  the  ambiguity  function  along  constant  range 
and  velocity  axes  are  linearly  related  to  the  group’s  first-  and  second-order  statistics  in  velocity  and 
position,  respectively. 


2 


WORK  COMPLETED 


The  research  work  in  Fiscal  Year  2008  was  a  great  success  in  developing  an  analytical  expression  for 
the  second-order  covariance  of  MLEs  based  on  general  complex  Gaussian  data  vectors  where  both  the 
data  mean  and  variance  are  parameter  dependent.  This  expression  was  then  used  to  calculate  necessary 
conditions  on  sample  size  and  SNR  for  accurate  source  localization  estimation  in  a  fluctuating  ocean 
waveguide.  We  also  quantified  the  effects  of  ocean  variability  by  calculating  and  comparing  the  biases 
and  errors  of  MLEs  of  source  position  both  in  a  static  waveguide  and  in  the  presence  of  internal  waves. 
We  investigated  the  advantages  of  using  the  derived  second-order  covariance  expression  by  repeating 
the  above  calculations  for  the  special  cases  where  either  the  covariance  or  the  mean  of  the 
measurement  are  assumed  to  be  parameter  independent. 

Also  in  Eiscal  Year  2008,  we  derived  an  analytical  expression  for  the  ambiguity  function  of  the 
acoustic  field  scattered  from  groups  of  underwater  targets.  This  expression  was  then  used  to  calculate 
the  moments  of  the  ambiguity  function  magnitude  squared  along  constant  range  and  velocity  axes, 
which  were  in  turn  shown  to  be  linearly  related  to  the  group’s  velocity  and  position  statistics,  given 
appropriate  signal  design.  Our  approach  fully  accounts  for  waveguide  dispersion  and  has  been 
developed  with  a  remote  long-range  platform  in  mind,  so  that  simultaneous  estimation  of  range  and 
velocity  over  continental  shelf-scales  may  be  possible.  Numerical  simulations  were  employed  to  model 
typical  cases  of  migrating  and/or  randomly  swarming  fish  schools  imaged  remotely  from  ranges  larger 
than  10  km. 

RESULTS 

We  developed  an  analytical  expression  for  the  second-order  covariance  of  MLEs  obtained  from 
general  complex  Gaussian  data  vectors.  We  used  this  expression  to  calculate  the  conditions  necessary 
to  obtain  source  localization  estimates  of  the  desired  accuracy  and  found  that  sample  sizes  and  SNRs 
may  increase  prohibitively  by  several  orders  of  magnitude  when  internal  wave  fluctuations  lead  to  a 
dominant  incoherent  intensity  component  in  the  total  acoustic  field  intensity.  We  also  investigated  the 
effect  of  assuming  that  the  measurement  corresponded  to  (i)  a  deterministic  signal  vector,  or  (ii)  a  fully 
randomized,  zero-mean  signal  vector,  both  embedded  in  additive  white  noise.  We  found  that  in  the  first 
case  the  CRLB,  as  well  as  the  sample  size  or  SNR  necessary  to  attain  it  are  significantly 
underestimated,  while  in  the  second  case  they  are  significantly  overestimated. 

We  developed  a  method  for  estimating  the  statistics  of  the  velocity  and  position  of  groups  of  self- 
propelled  underwater  targets  using  Ocean  Acoustics  Waveguide  Remote  Sensing  (OAWRS)  Doppler 
analysis.  We  demonstrated  the  method  using  illustrative  examples  for  simulated  shoaling  and 
migrating  fish  shoals,  as  well  as  for  a  group  of  AUV’s.  We  found  that  the  method  performs  very  well 
in  a  typical  continental  shelf  environment  and  that  it  is  often  possible  to  obtain  very  accurate  (i.e.  less 
than  10%  error)  estimates  of  the  velocity  and  position  statistics  for  a  single  target,  as  well  as  for  a 
group  of  targets.  We  also  found  that,  typically,  a  larger  number  of  targets  leads  to  more  reliable 
estimates.  Einally,  we  showed  that  the  second-order  statistics  can  provide  a  means  for  target 
discrimination  and  classification. 
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IMPACT/APPLICATIONS 


1,  Passive  Source  Localization  in  a  Fluctuating  Ocean  Waveguide 

The  calculated  biases  and  errors  may  be  used  to  quantify  the  effects  of  environmental  uncertainties  on 
passive  source  localization  techniques,  such  as  matched-field  processing  (MFP).  The  results  and 
methodology  presented  can  also  be  used  in  experimental  design  to  ensure  that  statistical  biases  and 
errors  meet  any  required  thresholds.  Further,  this  research  motivates  the  need  for  developing  analytical 
models  that  provide  analytical  expressions  for  the  statistics  of  the  acoustic  field  propagating  through 
random  and/or  fluctuating  environments. 

2.  Velocity  Estimation  of  Underwater  Target  Swarms 

It  was  recently  shown  that  densely  populated  fish  schools  may  be  the  dominant  cause  of  environmental 
clutter  in  continental  shelf  environments.^  Since  the  behavioral  dynamics  of  biological  scatterers  can 
be  expected  to  differ  from  those  of  man-made  targets,  it  is  worth  examining  whether  they  can  be 
inferred  and  then  used  as  a  means  of  clutter  identification.  Instantaneous  estimates  of  the  second-order 
statistics  of  both  the  velocity  and  position  of  large  groups  of  underwater  targets  may  then  provide  new 
opportunities  for  cutter  mitigation. 

PUBLICATIONS 
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•  Nicholas  C.  Makris  et  ah.  Critical  population  density  triggers  rapid  formation  of  vast  oceanic  fish 
shoals.  Science  323,  1734-1737  (2009). 

•  Nicholas  C.  Makris,  Purnima  Ratilal,  Michele  Zanolin  and  loannis  Bertsatos,  Obtaining  optimal  time- 
delay,  source  localization  and  tracking  estimates  in  free  space  and  in  an  ocean  waveguide,  J.  Acoust. 
Soc.  Am.  116,  2606  (2004). 

•  Nicholas  C.  Makris,  Yisan  Eai,  loannis  Bertsatos,  Deanelle  Symonds,  Sunwoong  Lee  and  Purnima 
Ratilal,  Expected  Doppler  shift  and  spread  in  long  range  scattering  from  fish  schools  on  the  continental 
shelf,  J.  Acoust.  Soc.  Am.  116,  2488  (2004). 

•  Purnima  Ratilal,  loannis  Bertsatos,  Tianrun  Chen,  Michele  Zanolin  and  Nicholas  C.  Makris,  Optimal 
passive  source  localization  in  a  fluctuating  ocean  waveguide  based  on  an  analytic  model  for  the  mean 
field  and  covariance,  J.  Acoust.  Soc.  Am.  115,  2473  (2004). 
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•  loannis  Bertsatos  and  Nicholas  C.  Makris,  Optimally  resolving  Lambertian  surface  orientation,  J. 
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